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Nonequilibrium  Superconducting  States 
with  Two  Coexisting  Energy  Gaps 


Gerd  Schbn^ 

Institut  ftir  Theorie  der  Kondensierten 
Mater ie,  Universitat  Karlsruhe,  BRD 


and 

Andre-M.  Tremblay 

Laboratory  of  Atomic  and  Solid  State 
Physics,  Cornell  University,  Ithaca 
New  York,  14853 


The  effect  of  tunnel  currents  in  superconducting 
junctions  on  the  energy  gap  is  calculated.  For  certain 
parameters  two  different  gaps  can  exist.  The  stability 
of  these  solutions  is  investigated  and  at  a  certain 

i 

voltage  a  "first  order  transition"  is  found.  This  result 
explains  the  experimentally  observed  inhomogeneous  states 
in  superconducting  tunnel  junctions. 


2 


Recent  experiments  on  superconducting  tunnel  junctions ^  ^ 

suggest  that  for  certain  injection  currents  and  voltages  a 

superconductor  sustains  simultaneously  two  different  values 

of  the  energy  gap.  Existing  phenomenological  models ^ 

predict  that  above  a  critical  density  of  the  excitations  nc, 

the  superconductor  has  an  intrinsic  instability  with  respect 

to  the  formation  of  a  spatially  inhomogeneous  state.  This 

is  not  in  agreement  with  the  experimental  results  of  Gray 
(2) 

and  Willemsen  .  These  authors  interpreted  the  effect 
by  inhomogeneities  in  the  probes,  with  a  lower  gap  value, 
which  grow  with  increasing  total  current. 

In  this  letter  we  describe  microscopically  a  super¬ 
conducting  tunnel  junction  consisting  of  an  injector  and  a 
probe  and  find  that  two  stable  values  of  the  energy  gap 
can  exist  in  the  probe.  In  accord  with  the  experiments, 
at  a  certain  voltage  a  first  order  phase  transition  takes 
place  where  the  part  of  the  probe  with  the  lower  gap  and 
larger  injection  current  density  grows  while  the  part  with 
the  larger  gap  and  lower  injection  current  density  decreases. 

The  relative  size  of  the  two  regions  is  controlled  by  the 
total  injection  current.  This  is  analogous  to  a  liquid- 
gas  transition  at  a  certain  pressure  where  the  relative 
volumes  are  controlled  by  the  total  volume.  Our  result  differs 
from  the  model  of  Ref .  2  insofar  as  we  find  the  existence 
of  two  gap  values  to  be  an  intrinsic  property  of  the  system. 
Furthermore,  our  approach  describes  the  gap  enhancement  by 

quasiparticle  extraction,  investigated  in  the  experiments 

( 4 ) 

of  Chi  and  Clarke  .  We  will  present  a  detailed  analysis 

p 

in  two  limits,  near  Tc  and  near  T  =  0.  Qualitatively,  the 
same  results  can  be  obtained  for  any  intermediate  temperature. 

We  assume  that  the  injector  is  thick  and  not  appreciably 
perturbed  by  the  current  and  that  the  phonons  remain  in 
equilibrium.  Thus,  the  probe  can  be  studied  using  the 
kinetic  equations  for  the  quasiparticle  distribution  function 
and  the  order  parameter  derived  in  Ref.  5.  In  general 


there  is  no  simple  relation  between  the  injected  current 
(branch  imbalance^)  and  the  total  number  of  excitations 
(electron- plus  hole- like)  created  in  the  probe.  Only  the 
latter  quantity  modifies  the  magnitude  of  the  order  parameter.  /  ■  >/ 
The  corresponding  deviation  of  the  distribution  function  5fg,  '  ■ 

which  is  odd  in  energy,  is  obtained  from  the  Boltzmann  / 

equation: 


M1(E)6£e  -  K(5f) 


PE  -  E)*3f°(E)/3E  •  AA/E 


v 

% 

■H,  % 


(E)  is  the  normalized  BCS  density  of  states  and  K(6f) 

Sv'./ 

describes  inelastic  electron-phonon  scattering.  It  suffices 
to  know  that  K(6f)  can  be  split  into  a  "scattering  out"  term 
-te  1  Jf£E)6£E,  and  a  "scattering  in"  term,  an  integral  operator. 
The  perturbation  is  PE  *  b  ^  (E)  J^V^fE  -  eV)(th  E/2T  -th(E  -  eV)/2T)  + 
-  eV)]  where  B”1  =  RJ28e2No,  R  is  the  resistance  of 
the  junction,  8  the  volume  of  the  probe,  N0  the  normal  density 
of  states, and  Af refers  to  the  injector.  The  thickness 
of  the  probe  is  small,  and  no  spatial  variations  in  this 
direction  occur.  For  the  moment  we  also  neglect  spatial 
variations  in  the  junction  plane.  The  effect  of  6f£  on 
A,  the  magnitude  of  the  order-parameter,  follows  for  tempera¬ 
tures  near  the  critical  temperature  of  the  probe  from  the 
Ginzburg-Landau  (G  -  L)  equation^ 


ffA/8T  -  xA 


(a  -  BA  ) A 


where  x  =  -  J dE  l/Ev^E)  6f£  and  a  *  (Tc-T)/Tc,B  =  7;  (3) /8it2T2 . 
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• 

We  obtain  analytic  results  if  we  assume  that  the  energy 
gap  of  the  injector  and  the  voltage  are  small  A^/T , | eV | /T  <<  1. 
The  linear  term  of  an  expansion  of  PE  in  powers  of  eV/T  is 
localized  in  a  narrow  energy  region  |e|  sAi  +  |eV|  <<  T. 


Consequently,  the  corresponding  contribution  6fE^  =  t0pe/<^i 

to  the  stationary  nonequi librium  distribution,  is  obtained 

neglecting  the  "scattering  in"  term  and  taking  the  scattering 

time  rE  at  E  *  0.  Depending  on  the  relative  values  of 

A, A.  and  V,  6f  J^describes  an  increase  or  decrease  of  the 
i  E 

density  of  excitations  near  the  gap  edge.  The  resulting 

X^(A,V)  also  changes  sign,  and  both  gap  reduction  or  gap 

enhancement  are  possible.  By  contrast,  the  quadratic  term 

of  the  expansion  of  P_  always  results  in  a  gap  reduction. 

(2)  * 

Since  6f_v  is  not  localized  near  the  gap  edge,  the 
£ 

"scattering  in"  term  cannot  be  neglected.  A  variational 
calculation^  yields  *  -1.4Btq  (eV/2T)2*  This 

contribution  can  be  interpreted  as  an  effective  increase  in 
temperature.  For  the  following  discussion  we  shall  neglect 
it,  although  it  may  be  important  for  a  quantitative  analysis. 

In  the  considered  limit  <A,A^, |eV|<<T) ,  we  find 

X(1)(A,V)  *  -20  (A  +  A±  -  |  eV  | )  B  T0|eV/2T(  g{  (d  -  |eV|)  K(k)  +  (c  -  d)  ir  (a2  ,k) 

(3) 

where  K  and  it  are  complete  elliptic  integrals  of  the  first 

and  third  kind,  g  *  2/  Ra  -  c)  (b  -  d)l  ,  a2  =  (b-c)/(b-d), 

2  2 

k  -a  (a-d)/(a-c),  and  a,  b,c,d  are  the  parameters  jeVl+A^nd 
±A  assigned  such  that  a>b>c>d.  In  order  to  obtain  the 
stationary  solution  of  the  G  -  L  equation  including  the  non¬ 
equilibrium  term,  we  employ  graphical  constructions.  The 

fl)  2 

intersections  of  ~x  (A,V)  and  the  curve  a  -  SA  yield 

possible  solutions  of  A.  A  large  positive  value  of  x^  (A,V) 

at  |eV|  *  A^  -  A  (for  A^ > A) ,  corresponding  to  a  net  extraction 

of  excitations,  results  in  a  large  gap  enhancement  at  this 

voltage.  This  is  in  qualitative  agreement  with  the  ex- 

14) 

perimental  results  of  Chi  and  Clarke  . 

For  our  present  problem,  larger  voltages  are  of  interest. 

As  shown  in  Fig.  1,  x^(A,V)  has  a  step  structure  at 
|eV|  »  A^  +  A.  This  step  corresponds  to  the  step  in  the 


T 


5  - 


I  (V)  characteristic  of  an  ideal  SIS  tunnel  junction  at  the 
same  voltage,  however,  is  finite  at  low  voltages 

|eV|  A+  A^  and  zero  above.  In  the  vicinity  of  the  step, 
we  can  approximate  Eq. (3)  by 

X(1)(A,V)  *  B  to  |eV/2T|ir/A^7TO<A+  A±  -  |eV|  )  .  (4) 

For  suitably  chosen  parameters  (e.g.,  eV  *  2.2A^  in  Fig.l) 
we  find  three  solutions  of  the  G  -  L  equation  denoted  by 
A^,  Aj/  A3*  (In  the  presence  of  any  level  broadening,  the 
step  actually  has  to  be  replaced  by  a  finite  slope,  therefore, 
A^  is  also  a  solution.)  In  addition  A  *  0  is  a  solution. 

From  an  analysis  of  the  time  dependent  equations,  we  find 
that  A  ■  0  and  A2  are  unstable,  whereas  A^  and  A3  are  locally 
stable.  At  low  voltages,  we  find  only  one  enhanced  super¬ 
conducting  solution,  whereas  at  high  voltages  we  find 
only  the  unperturbed  solution.  In  the  inset  of  Fig.l,  the 
solutions  A  as  a  function  of  the  voltage  are  shown. 

Obviously, in  the  range  where  A (V)  has  two  values  with 
Ai  4  Ai  >  |eV|  >Ai  +  A3,  two  significantly  differing  values 
of  the  current  density  are  obtained. 

In  order  to  find  which  of  the  two  locally  stable 
solutions  is  globally  stable,  we  follow  the  analysis 
performed  by  Schmid  '.  The  probability  of  a  solution  A  is 
given  by  W(A)a  exp  (-  ?"(A)  /T),  where  the  generalized  free 
energy  is  ? (A)  *  -2NQfl  ^  dA '  -  BA ' 2  +  x  (A ' 0  A '  .  is 

o 

plotted  in  the  inset  of  Fig.  1.  The  situation  is  clearly 
analogous  to  a  first  order  phase  transition.  At  low 
voltages,  A^  is  the  globally  stable  "phase".  At  a  certain 
voltage  VQ,  where  A^  +  A^>  |eVQ|  >  Ai  +  A3»  we  have 
£~(A^)  *  ?"(A3),  and  a  transition  between  A^  and  A3  can  take 
place.  With  increasing  total  current,  the  size  of  the 
region  in  the  A^  "phase"  with  small  injection  current 
density  decreases  in  favour  of  the  region  in  the  A3 
"phase"  with  large  injection  current  density.  At  high 


enough  voltages*  the  A3  "phase"  is  globally  stable.  In 

addition*  the  metastable  states  will  result  in  a  hysteresis 

(2) 

under  suitable  conditions.  The  effect  of  x  or  generally 
of  heating  will  be  to  reduce  A^  and  A^  below  the  values 
found  here*  but  these  processes  are  "smooth"  enough  that  the 
difference  A^  -  A^  will  be  preserved. 

All  of  these  results* including  the  dependence  of 
Ai ~  A,  on  the  junction  resistance,  are  in  good  agreement 
with  the  experiments.  Even  the  hysteresis  along  the  lower 
branch  of  I (V)  has  been  detected .  In  the  case  of  stimulation 
of  superconductivity  by  microwaves  or  by  phonons  at  a 

(8) 

certain  temperature  two  stable  solutions  are  found  too 
(normal  and  superconducting  or  both  superconducting  with 
different  gaps,  respectively).  However,  the  transition 
between  these  solutions  occurs  abruptly,  since  there  is  no 
external  variable,  such  as  the  total  injection  current  in 
the  present  problem  to  control  the  transition. 

Even  in  the  state  where  the  two  gaps  coexist,  the 
distribution  -function  <SfE  (in  the  limit  considered)  is 
single  valued.  Therefore,  spatially  inhomogeneous  problems 
are  described  by  adding  the  space  derivatives  to  the  G  -L 
equation.  We  find  ^  ^ that  at  V  *  VQ  a  stationary  wall  sep¬ 
arating  the  two  phases  can  exist.  Apart  from  the  fact  that  the 
location  of  the  wall  can  be  shifted  arbitrarily,  this 
solution  is  locally  stable.  This  confirms  our  result  of 
the  stable  coexistence  of  the  two  phases  at  the  voltage 
V  -  v  .  On  the  other  hand,  droplets  or  periodic  structures 
correspond  to  saddlepoint  solutions  and  are  therefore  un¬ 
stable. 

In  the  low  temperature  limit  (A,  A^  >>  T)  qualitatively 
similiar  results  can  be  found.  Particularly  interesting 
are  injection  voltages  close  to  the  sum  of  the  two  gaps: 

! (Aa  +  A)  -  )eV|  j  /  A  <<  1.  In  this  case  PE  and  consequently 
$f E  are  localized  near  the  gap  edge.  Integrating  the 
Boltzmann  equation  with  respect  to  the  energy,  we  find 
an  equation  for  the  total  number  of  excitations  per  unit 


volume  n  *  4Nq  ^dEA^CE)  5f£  in  the  form  of  a  linear  Rothwarf- 
0  (11) 

Taylor  equation'  .  In  the  stationary  case, 


2n/xR(E*A)  =  4Nq  ^  dEPE  *  jl(V)/en|  .  (5) 

o 

In  deriving  Eq.5,  the  localization  of  6f  allowed  us  to 

C* 

neglect  the  energy  dependence  of  the  recombination  rate  and 

take  its  value  at  the  gap  edge t  (E = A)  = 

—1  1/2  R 
6.5tq  (T/A)  '  exp(-A/T).  Furthermore,  within  terms  of 

order  exp  (-A/T) ,  J  dE  P£  is  proportional  to  the  injection 


current.  The  effect  of  nonequilibrium  excitations  on  the 
order  parameter  follows  from  the  self-consistency  relation 
A  »  XA  J^dEl/E  JfjiE)  (th  E/2T  -  26f£)  .  Again  the  localization 


o 

of  6fE  allows  us  to  simplify  this  equation  by  setting 
A/E  <5fE  c  6f£.  Thus,  to  lowest  order  in  e  =  (A-AQ)/A0, 
where  Aq  is  the  unperturbed  T  *  0  gap,  we  obtain 


c  =  -n/2NoAQ0((|eV|  -  Aq- Ai)/Ao- e)  .  (6) 

The  step  in  the  injection  current  at  the  voltage  |eV|  =*  A  +  A^ 
leads  to  a  corresponding  step  in  n,  which  is  made  explicit 
in  Eq. (6) .  For  suitably  chosen  parameters,  we  find  again 
two  locally  stable  solutions  e*0  and  e  =  -n/2NQA0.  The 
physical  interpretation  of  the  two  solutions  is  straight¬ 
forward.  First  we  notice  that  near  T  =  0  there  are  no 
excitations  which  could  be  extracted  from  the  system,  and 
an  injected  current  increases  their  number.  At  a  voltage 
|eV|  slightly  below  A^Aq,  the  probe  can  be  either  in  the 
unperturbed  gap  state  AQ  with  no  current  or  in  a  low  gap 
state  with  finite  injection  current  and  consequently  in¬ 
creased  excitation  number  stabilizing  the  low  gap  value. 


-  8  _ 


At  any  intermediate  temperature  o<T<T  ,  the  strong 
voltage  dependence  of  the  tunneling  current  at  |eV|  «  A^  + A 
results  in  a  steplike  modification  of  the  quasiparticle 
distribution.  For  this  reason  we  expect  the  qualitative 
result  of  this  paper,  the  coexistence  of  two  different  gap 
values,  to  be  true  at  all  temperatures  below  T^. 
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Figure  Caption 

Fig.  1:  Graphical  solution  of  the  G  -L  equation  for  voltages 

|eV|  >  A. .  The  constants  ct,B  and  B  in  are 

1  (2) 
chosen  arbitrarily,  the  higher  order  term  x  is 

neglected.  The  unperturbed  solution  is  chosen  to  be 

AQ  *  A^  X^(A,V)  has  a  step  at  A  =  |eV|  -A^  For 

the  same  parameters,  the  generalized  free  energy 

3^(A,VJ  is  shown  in  the  upper  inset,  while  the 

solutions  A(V)  are  shown  in  the  lower  inset. 


